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Abstract 

We study the dynamics of two-component Bose-Einstein condensates in periodic potentials in 
one dimension. Elliptic potentials which have the sinusoidal optical potential as a special case 
are considered. We construct exact nonstationary solutions to the mean-field equations of mo- 
tion. Among the solutions are two types of temporally-periodic solutions — in one type there 
are condensate oscillations between neighboring potential wells, while in the other the condensates 
oscillate from side to side within the wells. Our numerical studies of the stability of these solutions 
suggests the existence of one-parameter families of stable nonstationary solutions. 
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I. Introduction 

In conventional magnetic traps, the spins of the alkali atoms making up a Bose-Einstein 
condensate (BEC) are frozen, and all atoms are in the first hyperfine manifold. Thus, even 
if the atoms have spin, the condensate order parameter is a single complex scalar. 

Over the last few years, two methods to produce mixtures of two distinguishable BEC's 
have been developed and implemented. In these condensate mixtures, ip has two components. 
The Boulder group, for example, produced a mixture of two condensates consisting of two 
different hyperfine spin states of ^^Rb . Initially, a single condensate in the spin state 
|1 > with hyperfine spin F = 1 was trapped magnetically. A two-photon transition was 
then used to transfer a portion of the atoms to a second spin state |2> with F = 2. Once 
the two-photon pulse had ended, the number of atoms in each condensate was essentially 
constant. Modugno et al. on the other hand, produced a mixture of two Bose-Einstein 
condensates of different atomic species, ^^K and ^''Rb 0]. 

An atom placed in a standing light wave is subject to induced dipole forces. The resulting 
atomic potential is sinusoidal, and is referred to as an optical lattice. Merely by altering 
the phase, wavelength and intensity of the light, the position, lattice spacing and well depth 
of the optical lattice can be adjusted. Thus, optical lattices are well-characterized and 
controllable. 

By adiabatically increasing the depth of the lattice potential, trapped BEC's can be 
transferred into an optical lattice. BEC's on optical lattices have attracted a great deal of 
interest for several reasons. In experiments, the analogs of the ac and dc Josephson effects 

n n 

p], Bloch oscillations and Landau-Zener tunneling p, number-squeezed states j2|, and the 
Mott-insulator transition ^ have been studied. In addition, an optical lattice in which each 
site is occupied by a single alkali atom in its ground state is a promising candidate for a 
register in a quantum computer. One route currently being explored to initialize such a 
register is to begin with a BEC and then adiabatically turn on an optical lattice, with the 
result that one atom in its ground state occupies each site From a theoretical standpoint, 
BEC's on optical lattices are interesting because the simple sinusoidal form of the potential 
has made it possible to obtain some exact solutions to the mean-field equations of motion, 
in spite of their nonlinearity jiol . 11, 12 , 13 1 . 



Two condensates confined to a cigar-shaped trap are quasi-one-dimensional if the trap is 
highly elongated and if the tranverse dimensions of the trap are comparable to the healing 
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lengths. A quasi-one-dimensional optical lattice can be realized by superimposing a standing 
light wave on the cigar-shaped trap. In the mean-field approximation, the condensates are 
described by coupled nonlinear Schrodinger equations in one space dimension with periodic 
external potentials. 

In this paper, we study the dynamics of two-component condensates in periodic potentials 
in one dimension (ID). We will consider elliptic potentials which have the sinusoidal potential 
as a special case. Exact stationary solutions to this problem have been found by Deconinck 
et al. We will use these solutions to construct exact nonstationary solutions to the 

problem. Our numerical studies of the stability of these solutions strongly suggest that there 
are one-parameter families of stable solutions. 

We will only study condensate mixtures in which the atoms of the different components 
have the same atomic mass /x. This means that mixtures of condensates of different atomic 
species are excluded from our analysis, and accordingly we will confine ourselves to consid- 
eration of the two-component condensates like those studied by the Boulder group 

The paper is organized as follows. In Section II, we introduce the mean-field equations of 
motion and develop a general method of constructing nonstationary solutions. This method 
is applied to condensates in a ID elliptic potential in Sec. Ill, and the physical interpretation 
of the resulting solutions is discussed in Sec. IV. The stability of the solutions is the subject 
of Sec. V. Finally, we summarize our results in Sec. VI. 

II. General Development 

In the two-component condensates studied by the Boulder group, the number of atoms 
in each condensate is nearly constant jjj, |2| . We will therefore neglect loss of atoms to the 
normal phases and switching of atoms from one condensate to another. The mean-field 
equations of motion are then 

^^^ = -^^+(E«.IV^d^j^. + ^.V^. (1) 

for j = 1, 2. Here ipj = tpj{x,t) is the condensate wave function for the jth species and aji 
describes the interaction of an atom in the jth condensate and an atom in the Zth condensate. 
The 2x2 symmetric matrix M = {aij} will be referred to as the interaction matrix. The 
interaction strengths an, aio and 022 are known to the 1% level for ^'^Rb, and are in the 
proportion 1.02 : 1 : 0.97 We shall neglect the difference between the interaction 

strengths and set a = an = ai2 = a22- 



In the case of an optical lattice, the potentials Vi and V2 are sinusoidal. Vi and V2 are 
very nearly equal for linearly polarized light, provided that the detuning is not too small 
We shall neglect the difference Vi — V2 and set 

V = Vi = V2 = -Vo sin^ qx, (2) 
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where Vq is a constant which depends on the intensity of the light, q = 27r/A and A is the 
optical wavelength. In fact, we will study the more general potential 

V = -Vosn\qx,k), (3) 

where sn {qx, k) is the Jacobian elliptic sine function with elliptic modulus k G [0, 1]. The 
potential (jH)) reduces to the optical potential Q for = and to the single potential well 
(or barrier) V = —Vq tanh^ qx for k = 1. For < /c < 1, the potential is periodic with period 
2K{k)/q, where 

-'0 VI - k^ sm^ 9 |— I 

Plots of the potential Q for a range of k values may be found in Ref. To simplify the 

notation, we will not display the k dependence of the Jacobian elliptic functions from this 
point on. 

Let ip = {1^1,1^2)'^ ■ With the simplifications we have made, the equation of motion is 

= + Vosn^iqx) ip. (5) 

This equation is further simplified by introducing the dimensionless variables x = qx, t = 



hqH/fi, ipj = {^\a\fi/hq)ipj, a = sgn(a) and Vq = fiVo/{hq)^ and then dropping the tildes. 
This gives 

dip 1 d'^ij) 

i-g^ = -^-g^ + aiip^i/j)ip -Vosn'^{x)ip. (6) 



If f/ is a 2 X 2 unitary matrix, 

if' = Uip (7) 

is also a solution to the equation of motion (jH} This observation will allow us to 

construct nonstationary solutions from the stationary solutions of Deconinck et al. |12l • 
III. Construction of the Nonstationary Solutions 

Nonstationary solutions can be constructed for both a = +1 and a = —1. For a = —1, 
any two atoms attract each other, regardless of whether or not they belong to the same 



condensate. Thus, the nonstationary solutions are Ukely unstable against collapse for a = 
— 1. We will therefore restrict our attention to the case a = +1 for the remainder of the 
paper. 

Let A and B be arbitrary nonnegative real numbers. From Ref. 
types of stationary solutions to Eq. (jH)): 
Type I with 



12 1 we obtain three 



A cnx 



exp 



and 



Type II with 



t{- + A^ + B'']t 



(8) 



(9) 



/ Asnx \ 



[i + e) + A^ + B^' 



and 



V={k'-A' + k'B')sn'x + A 



(10) 



(11) 



and lastly Type III with 



Asnx 



kB cnx e*'^^*/^ 



exp < —I 



[I + F) + A' + k^B^ 



(12) 



and 



V={e-A^ + k^B^)sn^x + A^. (13) 

We now apply a unitary transformation to these solutions: we set ip' = U{6)iIj, where 

^ cos 6 sin 6 \ 
V — sin cos 6 I 



U{9) 



(14) 



Formally, we can think of ip as being the wave function of a spin-1/2 particle. U{6) then 
represents a rotation in spin space through an angle f3 = —26 about the y-axis. The external 
potential is unchanged by the transformation (l8l |. 

We may confine our attention to the interval < < tt because U{tt)iP' = In 
fact, it is possible to reduce this interval. To see this, let ^{x,t;9) = 'ip'{x,t) and note 
that ^(x, t;e + 7r/2) = [/(7r/2)^(x, t; 9) = {^2{x, t; 6), -^i{x, t; e)f. Now if (V^i, V^a)^ is a 
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solution to the equation of motion (jH)), then so is (e^'^^Vi) e^'^'^ip2Y' for arbitrary real constants 
(pi and 02- Moreover, these solutions have the same physical meaning. Apart from a switch 
in the labels 1 and 2 and an irrelevant phase change, therefore, ^(x, t; + 7r/2) and ^(x, t] 9) 
are identical, and it suffices to consider the interval < 9 < -k /2. There is a symmetry 
about 9 = Ti/A that enables us to further pare down the range of 9 we must consider. We 
begin with the observation that 

For Types I, II and III, ip has the form (/i(a;)e~*'^^*, /2(x)e~*'^2*)^. Employing the identity 
fITSll. we obtain 



^ ( X, t; ^ - ) = exp 



UJ2 

-ITX ' 



° ^]^(x,t ^-,9]. (16) 

10/ V UJ1-U2 



.UJ2 - tOl- 

This means that except for a switch in the labels 1 and 2 and an irrelevant phase factor, 
\E'(x, t]^ — 9) and \E'(x, t — ; 9i) are the same. Thus, we need only consider the interval 
< < 7r/4. 

= U{9)ip is a nonstationary solution for < < 7r/4 if t/^ is a stationary solution of 



Type I, II or III. Other types of stationary solution are constructed in Ref. Il2|: There are 
solutions in which both ipi and ip2 are proportional to the same Jacobian elliptic function, 
as well as solutions in which is a linear function of a Jacobian elliptic function for i = 1 
and 2. However, when a unitary transformation is applied to solutions of these types, the 
result is simply another stationary solution of the same kind. Therefore, we will not consider 
these solutions further. 

IV. Physical Interpretation of the Nonstationary Solutions 

A. The Trigonometric Limit k = 

In the k = limit, the Jacobian elliptic functions reduce to trigonometric functions. For 
solutions of Type I, the density of condensate 1 is 

n[ = = cos'^ 9 cos^ X + sin"^ 9 + AB sin 29 cos xcos{t/ 2), (17) 

while the density of condensate 2 is given by 

n'2 = |V^2p = A^sin^^cos^x + E^cos^^- A5sin2^cosxcos(t/2). (18) 
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n'^ and n'2 oscillate in time with period T = Air (in unsealed physical units T = A-Kfi/hq"^). 
The potential V = sin^ x is an optical potential with minima at the points x = mir, 
where m is any integer. We divide the lattice of potential minima into two sublattices: 
sublattice 1 with even m, and sublattice 2 with odd m. The total condensate density 
n' = n[ + ^2 = cos^ X + B"^ is independent of time, and its maxima occur at the potential 
minima. At time t = 0, the global maxima of n'l are on sublattice 1, while at time t = T/2, 
they are on sublattice 2 (Fig. 1). Naturally, at time t = T, the global maxima of n[ have 
returned to sublattice 1. Note that the global maxima reside on one of the two sublattices 
for all times t, and their locations change discontinuously as time passes. The global maxima 
of ^Iso oscillate between sublattices 1 and 2, but the oscillations of lag those of n[ by 
half a period. When the condensates move between the two sublattices in this way, we will 
say that they undergo sublattice oscillations. 

In the k = limit, the solutions of Types I and II become identical, apart from a 
translation of both ip' and V through 11/2. The Type III solution reduces to a stationary 
solution already studied by Deconinck et al. 3] and will not be considered further here. 

B. Solutions with < A; < 1 

We now turn to the nature of the solutions for < A; < 1. Explicitly, the Type I solution 

is 2 

/ A cos6'cnx + -B sin6'dnxe*'^^~'^ \ r /I \ 1 

^'=( . 2^ / exp -z + + . (19) 

V — v4 sin 6'cna; + -B cos 6'dna; e*^^ ^1^1'^) L V2 /J 

The potential minima appear on the lattice of points x = 2'mK, where m is any integer. 
We again divide the lattice of potential minima into two sublattices: sublattice 1 with even 
m, and sublattice 2 with odd m. The two condensates execute sublattice oscillations with 
temporal period T = An/ {1 — k"^), as shown in Fig. 2. 

The nature of the solutions of Types II and III is more complex. We will continue to 
refer to the set of points x = 2mK with integer m as the lattice, and to divide this lattice 



into sublattices 1 and 2. The potential minima are on the lattice for < A < kVl + W. 
For A > k^\ + on the other hand, the potential maxima are on the lattice, and the 
minima lie on the set of points x = IK, where I is any odd integer. 

Let < n'j > be the temporal average of n'j. For solutions of Types II and III, we say that 
the motion of condensate j is of type a (/5) if the maxima of < n'j > are located at the odd 
(even) multiples of K. The motion of the two condensates will be said to be of type 7^ if 
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the motion of condensate 1 is of type 7 and the motion of condensate 2 is of type 5, where 
7 and 5 can be either a or /3. 

We begin by considering the solutions of Type II. If the motion of condensate j is of type 
f3, the maxima of n'j smoothly oscillate from side to side as time passes and the oscillations 
of adjacent maxima in n'^ are 180° out of phase. Let Ni''^ be the number of atoms of species 
j on the interval {I — 1)K < x < {1 + 1)K for all odd integers 1. If the motion of condensate 
j is of type a, then N^2.+i oscillates periodically in time for each integer m. N^2.-i ^^^o 
oscillates periodically in time, but the oscillations of N^2.-i ^^S those of A^i^+i by half a 
period. For motion of both types a and the period of the temporal oscillations T is An. 

As shown in Fig. 3, motion of types aa, af3, and f3f3 occurs for different ranges of 
the parameters 9 and A/{kB). (Motion of the type (3a occurs in part of the region with 
7r/4 < 9 < 7r/2; these solutions are mapped to the a (3 sector in the region with < 9 < tt/A 
by the symmetry transformation discussed in Sec. III.) In the sector of the "phase 
diagram," the maxima of both n[ and ^2 reside in the potential wells. The maxima of n[ 
oscillate from side to side within the wells, and the oscillations of adjacent maxima are 180° 
out of phase (see Fig. 4). The maxima of 77,2 oscillate in the same fashion, but the temporal 
oscillations of 77,2 lag those of n'l by half a period. 

If A/{kB) is greater than both cot^ and Vl + the motion is of type aa, and the 
potential minima are at the points x — IK, where I is odd. The number of atoms of 
condensate 1 in the potential well centered on the point x — (4m + 1)K (i. e., A^im+i) 
oscillates periodically in time for each integer m. The number of atoms of condensate 1 
in the potential well centered on the point x = (4m — l)K (i. e., N^il^i) also oscillates 
periodically in time, but the oscillations of A^i^-i lag those of A^im+i by half a period. The 
motion of condensate 2 is analogous to that of condensate 1, except that it lags that of 
condensate 1 by half a period. Finally, the total condensate density n' — n[ + n'2 depends 
on position but is independent of time. 

For k <1/ -\/2, the motion just described is simply a sublattice oscillation, and the global 
maxima of n[ and n'g are always located at potential minima. A curious but interesting 
type of motion can occur for A; > 1/ \/2, though. In this case, for certain values of the 
parameters, as a global maximum in n'j grows in amplitude, it is initially located at the 
potential minimum. However, as its height continues to increase, the global maximum can 
split into two global maxima that are located to either side of the potential minimum (see 
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Fig. 5). 

In the region of the phase diagram in which 1 < A/{kB) < and < 6 < tt /4, 

motion of types aa and a/5 occurs, but the maxima of the total condensate density n' are 
located at maxima of the potential. These solutions are therefore expected to be unstable, 
and this is confirmed by our numerical simulations. 

We now turn to the nature of the solutions of Type III. In motion of both types a and 
f5, the maxima and minima of n'j oscillate continuously from side to side, and all of maxima 
and minima remain in phase with one another. The phase diagram for the solutions of 
Type III is given by Fig. 3, just as it was for the solutions of Type II. In the region of the 



phase diagram in which 1 < A/{kB) < ^/l + B^ and < 6 < n/A, the maxima of the total 
condensate density n' are located at maxima of the potential, and the solutions are expected 
to be unstable 19|. Outside this region, the maxima of n' coincide with the minima of the 
potential. If A/{kB) is greater than both cot 6 and Vl + B~'^, the motion is of the type aa 
and the maxima of both n'l and n'g reside in the potential wells. The maxima of n'l oscillate 
continuously from side to side and in phase with one another (see Fig. 6). The oscillations 
of the two condensates are 180° out of phase, and, as a result, the total condensate density 
n' does not depend on time. Qualitatively speaking, the same type of intra-well oscillation 
occurs in the f3(3 sector of the phase diagram. The two condensates oscillate in phase with 
one another in the ajS sector of the phase diagram, but the maxima of < n[ > coincide with 
the minima of < n'2 >. For all solutions of Type III, the period of oscillation T = Air/k'^. 

For solutions of both Types II and III, an interesting special case is obtained for A/(kB) = 
VT+~B~^. The external potential V is simply a constant in this case. Motion of type a(3 



occurs if 6* < cot^^(vT+5^); otherwise, the motion is of type aa. 
C. The Hyperbolic Limit k = 1 

The Jacobian elliptic functions reduce to hyperbolic functions for = 1. The Type I 



solution becomes a stationary solution studied in Ref. [1^. The Type II and III solutions, 
on the other hand, are identical for k = 1. Since these solutions are nonstationary, we will 
briefly consider their physical interpretation. 

The external potential for Type II solutions is ^ = {1 — + i?^)tanh^a; + A"^. A 
method of producing a potential with a tanh^ x spatial dependence in an experiment is yet 
to be found. The case A = vT+~B^ jg therefore of particular interest, since V is simply 
a constant in that case. For cot-i(Vl + 5-2) < 9 < 7r/4, a bound pair of dark solitons 
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oscillate about the origin. The two dark solitons move 180° out of phase with one another. 



For < < cot^^(Vl + B~'^), on the other hand, a dark soliton and a bright soliton oscillate 
in phase about the origin. The bright soliton cannot exist in isolation since the atoms repel 
one another, but the presence of the dark soliton stabilizes the bright soliton. Both of these 
kinds of solution — oscillatin g da rk-dark and dark-bright soliton pairs — have previously 
been found by Park and Shin |lq . 

This discussion suggests an alternative way of thinking about the special cases mentioned 
at the close of Section IV C. For both Types II and III, solutions of type a/5 have dark-bright 
soliton pairs oscillating about the lattice points, while solutions of type aa have bound pairs 
of dark solitons oscillating about these points. The dark-bright pairs oscillate in phase with 
each another. In contrast, the oscillations of the pairs of dark solitons are 180° out of phase. 
Neighboring solitons in n'j oscillate 180° out of phase in Type II solutions, whereas for Type 
III solutions neighboring solitons in n'^ oscillate in phase with one another. 

V. Numerical Investigation of the Stability of the Nonstationary Solutions 

tjj' = U{9)ip is a stable solution to the equation of motion (jHl) if and only if ^/S* is a 
stable solution to that equation. This observation has two notable consequences. First, 
analytical results on the stability of stationary solutions were obtained in Ref. and if is 
a stationary solution. Unfortunately, these analytical results are of no use here, as they only 
apply if the interaction matrix M is nonsingular. We will therefore probe the stability of 
the nonstationary solutions numerically. Secondly, in our numerical work it is sufficient to 
investigate the stability of our solutions for a single value of 6, which we choose to be ir/A. 

Our numerical method consists of solving the equation of motion ^ using a fourth-order 
Runge-Kutta method in time with a filtered pseudo-spectral method in space. For each run, 
a small amount of white noise was added to the initial data and the time for the onset of 
instability t* was determined. For all simulations, 

. /, \Un{t)-Uit)\ \ 

t = mm < t : — -— = 0.1 > , 

*>0 I Unit) + U{t) J 

where U{t) = max^{\ilj[{x,t)\,\ip'2{x,t)\}, f/„(t) = max^ {\'ip[^{x,t)\,\'ip2nix,t)\}, and ip'^^ 
and 1^271 represent the first and second components of the numerical solutions for and 1^2, 
respectively. 

The parameter space of our nonstationary solutions is too large to permit a comprehen- 
sive study of all possibilities. Instead, we will consider specific one-parameter families of 
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solutions, and determine how t* varies as the parameter changes. 

A. Trigonometric solutions close to the uniform solution 

The class of solutions considered first is trigonometric, and hence k = 0. Setting 6 = 7r/4, 
A = ^/2e, B = 1/2 and changing the zero of energy, the solution of Type 1 becomes 

/I + ee"**/^cosa;\ 

^'=1 -u/2 (20) 
\ 1 — ee ' cosx / 

with 

y = -2 - 2ecos2x. (21) 

As the parameter e — 0, this solution approaches the stationary uniform solution, which is 

1 i 

known to be stable |12]. We investigated the stability of this class of solutions for a range of 
values of e. In Fig. 7, t* is plotted as a function of e. The instability onset time t* appears 
to diverge as e approaches e'^ from above, where ~ 0.2. Two numerical runs are not 
displayed in Fig. 7: we find that t*{e = 0.205) > 60000 and t*{e = 0.2) > 132000. These 
results strongly suggest that t* is infinite for e < e^, and hence that there is a one-parameter 
family of stable solutions given by Eq. (pUj) with e < e'^ ~ 0.2. 

The behavior of an unstable solution with e = 1/2 is illustrated in Fig. 8. For this value of 
e, the solution is far from being uniform even at early times. After the onset of the instability 
at time t* ^ 450, the solution loses much of its coherence, although some is retained and 
new, larger-scale structures develop, as seen in Figs. 8 (c) and (d). 

B. Trigonometric solutions with increasing offset 

The second class of solutions considered is also trigonometric. With 6 = tt/A, A = 1, 



B = J —\ — 2e and V = 2t + \ — cos^ x, the solution of Type I becomes 



1 , e-'*/2 s 

, „ -I - e H 7^ cosx\ 

In general, we define the offset of a solution to be the smallest value that min[n'^(x, t), n2(x, t)] 
takes on for all x and t. The offset of the solution (1221) increases as e ^ — 00. It was shown 

n 

in Ref. [12] that increased offset stabilizes some stationary solutions. This motivated us 
to consider the class of solutions (j^^ with increasingly more negative e. In Fig. 9, the 
instability onset time t* is shown for a range of e values, t* remains finite for — e as large 
as 10, and it seems likely that t* is finite for all — e < 00. If this is indeed the case, the 
solutions (1^ are all unstable. 



11 



C. A class of elliptic solutions 

The third class of solutions we consider is not trigonometric, i. e., 7^ 0. In all numerical 
runs, k = 0.999, and so the solutions and the external potential are far from trigonometric. 
With 6 = 7t/A, A = e and B = 1 + e/k, the solution of Type II is 



i)' =\ ^ V N , (23) 

V -^e-*^i*snx + -ij (1 + I j e-*^2tdnx / 

with tui = (1 + A;^)/2 + (1 + e/A;)^, ^012 = A;^/2 + (1 + e/A;)^, and an appropriate additive 
shift of the potential. As e ^ 0, this solution approaches a stationary solution which has 
been proven to be linearly stable in the case of a nonsingular interaction matrix This 



motivates the consideration of the family of solutions ()23|) with parameter e. In Fig. 10, 
the instability onset time t* is shown for different values of e. The behavior of t* is quite 
interesting, as it displays one local maximum and one local minimum in the range of e 
considered. It is possible other local extrema exist for values of e > 1.5, but solving the 
equation of motion numerically becomes progressively more difficult as e increases because 
the nonlinearity grows stronger. As e approaches e'^ ~ 0.375 from above, t* appears to 
diverge. Thus, there seems to be a whole band of stable nonstationary solutions with the 
stationary solution clS cl limiting case. 

For values of e > e'^, the solution is unstable, and interesting structures appear 
after the onset of the instability. Figs. 11 - 13 illustrate various aspects of this behavior 
for k = 0.999 and e = 1/2. From these figures, we see that the solution becomes unstable 
shortly before t = 1640, and then is modulated in time with a period long compared to that 
of the intra-well oscillations. The amplitude of the density oscillations within the potential 
wells varies by up to a factor of five. After about three periods of the modulation, the 
modulated state in turn becomes unstable shortly after time t = 2600. 

These phenomena are also readily observed in the Fourier transform of the solution, which 
is shown in Fig. 13. Note the change in the spectrum during the modulated stage of the time 
evolution: the amplitudes of the Fourier coefficients change, but the number of activated 
modes remains approximately constant. The modulated state itself becomes unstable after 
a few periods of the modulation, and a band of approximately 60 unstable modes is excited, 
resulting in an apparent loss of coherence in the density. 
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VI. Summary 

In this paper, we studied the dynamics of two-component Bose-Einstein condensates in 
periodic potentials in one dimension. EUiptic potentials which have the sinusoidal optical po- 
tential as a special case were considered. We constructed exact nonstationary solutions to the 
mean-field equations of motion by performing a unitary transformation on previously-known 
stationary solutions. Among the solutions are two types of temporally-periodic solutions — 
in one type of solution there are condensate oscillations between neighboring potential wells 
and in the other type the condensates oscillate from side to side within the wells. Our 
numerical studies of the stability of these solutions suggests the existence of one-parameter 
families of stable solutions for both sinusoidal optical potentials and for elliptic potentials. 

R. M. B. thanks J. Roberts for useful discussions; B. D. and J. N. K. acknowledge support 
from the National Science Foundation (DMS-0139093 and DMS-0092682 respectively). 
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Figure Captions 

1. n\ (shown in dark gray) as a function of space and time in the optical potential (shown 
in light gray). The parameter values for this case are k = A = B = l/\/2 and 9 = tt/8. 
Note that the temporal oscillations of n'g lag those of n[ by half a period. 

2. n[ for a solution of Type I as a function of space and time. The external potential 
V is shown at the rear of the figure. The parameter values for this case are k = 0.99, 
A = B = I/V2 and 9 = n/8. As before, the temporal oscillations of n'g lag those of n[ by 
half a period. 

3. A "phase diagram" showing the types of motion which occur for different values of 
the parameters 9 and A/kB. This figure applies to solutions of both Types II and III. For 
< A/kB < Vl + 5-2, the potential minima are on the lattice of points x = 2mK, where 
m is an integer. For A/kB > Vl + on the other hand, the potential maxima are on 
the lattice, and the minima lie on the set of points x = IK, where / is any odd integer. 

4. n'l for a solution of Type II in the f3f3 sector of the phase diagram as a function of 
space and time. The external potential V is shown at the rear of the figure. The parameter 
values for this case are k = 0.99, A = 1, B = 3 and 9 = 71/ 4. Note the oscillations of the 
condensate from side to side within the potential wells. 

5. n[ for a solution of Type II in the aa sector of the phase diagram as a function of 
space and time. The external potential V is shown at the rear of the figure. The parameter 
values for this case are k = 0.99, A = 2.5, B = 1 and 9 = it /A. Note the periodic splitting of 
the peaks in n[ as the overall amount of condensate 1 in the potential wells rises and falls. 

6. n[ for a solution of Type III in the aa sector of the phase diagram as a function of 
space and time. The external potential V is shown at the rear of the figure. The parameter 
values for this case are k = 0.99, A = 3, B = 1 and 9 = vr/S. 

7. The instability onset time t* for the solutions given by Eq. ()20p. t* appears to diverge 
as e approaches ~ 0.2 from above. 

8. The onset of the instability for the trigonometric solution (|20|) with e = 1/2 in a 
computational domain of size L = Svr. Figures (a) and (b) show the density of the first and 
second condensate components, respectively. Figures (c) and (d) are the corresponding gray 
scale plots. Regions of high (low) density are shaded black (white). 

9. The instability onset time t* for the solutions given by Eq. ()22|) . 

10. The instability onset time t* for the solutions given by Eq. (|23|) . t* appears to diverge 



15 



as e approaches ~ 0.375 from above. 

11. The behavior of the elliptic solution ^ with e = 1/2, A; = 0.999, L = 167^(0.999) 
and ^ = X + L/16, immediately after the onset of the instability. Figures (a) and (b) show 
the densities of the first and second condensates, respectively. Figures (c) and (d) are the 
corresponding gray scale plots of these quantities. 

12. The behavior of the elliptic solution ^ with e = 1/2, k = 0.999, L = 16K(0.999) 
and C, = X + L/16, from before the onset of the instability until well after the end of the 
modulated phase of the motion. The left and right figures are gray scale plots of the densities 
of the first and second condensates, respectively. 

13. Numerical Fourier spectrum of the solution (|23p with e = 1/2 and k = 0.999. Only 
60 of the 512 modes used are shown. The left and right figures are the Fourier transforms 
of n[ and n'j, respectively. 
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